Cosmological Evolution of Interacting Phantom Energy with Dark Matter 
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We investigate the cosmological evolution of an interacting phantom energy model in which the 
phantom field has interaction with the dark matter. We discuss the existence and stability of scaling 
solutions for two types of specific interactions. One is motivated by the conformal transformation 
in string theory and the other is motivated by analogy with dissipation. In the former case, there 
exist no scaling solutions. In the latter case, there exist stable scaling solutions, which may give a 
phenomenological solution of the coincidence problem. Furthermore, the universe either accelerates 
forever or ends with a singularity, which is determined by not only the model parameters but also 
the initial velocity of the phantom field. 

PACS numbers: 98.80.Cq, 98.80.-k 



I. INTRODUCTION 

Scalar field plays an important role in modern cosmol- 
ogy. The dark energy can be attributed to the dynamics 
of a scalar field, for instance quintessence 0, Q, which 
convincingly realize the present accelerated expansion by 
using late-time attractor solutions, in which the scalar 
field mimics the perfect fluid in a wide range of parame- 
ters. Much attention has been drawn to the case of expo- 
nential potentials. The exponential potentials allow the 
possible existence of scaling solutions in which the scalar 
field energy density tracks that of the perfect fluid (so 
that at late times neither component can be negligible). 
In particular, a phase-plane analysis of the spatially flat 
FRW models showed that these solutions are the unique 
late-time attractors whenever they exist 0, More- 
over, exponential potentials appear naturally in the low 
dimensional models of string/M-theory 

The recent SNe data seem to favor the dark energy 
with the present equation of state w < — 1 |(| . To obtain 
w < — 1, phantom field with a negative kinetic term may 
be a simplest implementing and can be regarded as one of 
interesting possibilities describing dark energy [?]. The 
physical background for phantom type of matter with 
strongly negative pressure would be found in string the- 
ory [8j- Phantom field may also arise in higher-order 
theories of gravity 0, Brans-Dicke and non- minimally 
coupled scalar field theories . The cosmological mod- 
els which allow for phantom matter appear naturally in 
k-essence models In spite of the fact that the field 
theory of phantom fields encounters the problem of sta- 
bility which one could try to bypass by assuming them 
to be effective fields it is nevertheless interesting 

to study their cosmological implication. Recently, there 
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have been many relevant studies on this topic |13| . 

The physical properties of phantom energy are rather 
weird, as they include violation of the dominant energy 
condition and increasing energy density with the expan- 
sion of the universe. The latter ultimately leads to an un- 
wanted future singularity called big rip. This singularity 
is characterized by the divergence of the scale factor in 
a finite time in future [T^j . To avoid the cosmic dooms- 
day, some phantom field models were proposed ^ 
requires a special class of phantom field potentials with 
a local maximum. Moreover, the energy density of the 
phantom field increases with time, while the energy den- 
sity of the matter fluid decreases as the universe expands. 
Why are the energy density of dark matter and the phan- 
tom energy density of the same order just at the present 
epoch? This coincidence problem becomes more diffi- 
cult to solve in the phantom model. Therefore, from this 
point of view the cosmological scaling solution would be 
desirable for the history of the universe. Throughout this 
paper we use "scaling solution" as a meaning that the en- 
ergy densities of the phantom field and the dark matter 
are proportional. However, as shown in Ref. [16j, there 
exist no scaling solutions because the phantom energy 
increases while the matter energy decreases with time. 
But in the presence of the suitable interaction this case 
could be realized easily in Ref. 0. In this paper we in- 
vestigate the stability and existence of scaling solutions 
in the scenario of interacting phantom energy with dark 
matter. We consider two phenomenological models. One 
is motivated by the conformal transformation from the 
Jordan frame to Einstein frame in string theory and the 
other is motivated by analogy with dissipation. In the 
former model there exist no scaling solutions. However, 
in the latter model a phase-plane analysis shows that 
there exist two kinds of stable scaling solutions, which 
lead to two different fates of the universe. The universe 
either accelerates forever or ends with a big rip, which is 
determined by not only the model parameters but also 
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the initial velocity of the phantom field. We will also 
discuss the physical consequences of these results. 

II. INTERACTING PHANTOM ENERGY WITH 
DARK MATTER 



where the prime denotes a derivative with respect to the 
logarithm of the scale factor, N = In a, and the Fride- 
mann constraint equation Q becomes 

-x 2 +y 2 + z 2 = 1. (11) 



Let us consider a universe model where both the phan- 
tom field <p an d the dark matter p m are present. The 
Fricdmann equation in a spatially flat FRW metric can 
be written as 



The critical points, where x' = 0, y' = and z' = 0, 
correspond to an expanding universe with a scale factor 
a(t) given by a tx t 2 ^ 3 ( 1+WeS » . The effective equation of 
state for the total comic fluid is 



H = y(Pp + Pm), 



(1) 



where k 2 = SttGn is the gravitational coupling and the 
energy density and pressure, p p and P p , of the homoge- 
neous phantom field (f> are given by 

1 



Pv = "2< 

Pr, = 



(2) 
(3) 



respectively, in which V(4>) is the phantom field poten- 
tial. We postulate that the two components, p p and p m , 
interact through the interaction term Q according to 



p m + 3H(p m + P m ) = Q, 
p p + 3H{p p + P p ) = -C 



(4) 
(5) 



Suppose the dark matter possesses the equation of 
state P m = 0. The dynamics of the phantom field with 
an exponential potential 

7(0) = V exp(-A^) (6) 

has been analysed in Ref. • We assume the dimension- 
less constant A is positive since we can make it positive 
through 4> — > — if A < 0. Here we generalize the analy- 
sis to the case in which the phantom field has interaction 
with the dark matter. As we will see, this gives rise to 
some interesting novel features. Following Ref. 0, we 
define the following dimensionless variables 



K(j) 

Tiff' 



v 



V3H 



V3H 



(7) 



Notice that x 2 , y 2 and z 2 give the fraction of total en- 
ergy density carried by the field kinetic energy, the field 
potential energy and the dark matter, respectively. Thus 
the fractional densities of p p and p m can be written as 
il p = —x 2 + y 2 and fi m ~ z 2 , respectively. The evolution 
equations and JSJ can be written as the following set 
of equations: 




(12) 



Interaction terms Q have been discussed in the liter- 
ature within the context of inflation and reheating. In 
the conventional reheating model, an interaction term 
r^ip 2 dominates at the end of inflation when the scalar 
field is oscillating about the minimum of its potential. 
During this reheating phase the energy transferred from 
the scalar field is completely converted into the matter. 
Within the context of exponential potentials, an interac- 
tion term of the form Q = —en p m 4> was considered in 
Ref. [3[l3. It was shown that the matter scaling solu- 
tions were stable solutions and the age of the universe is 
older when such the interaction term is included. Certain 
string theories in which the energy sources are separately 
conserved in the Jordan frame naturally lead to interac- 
tion terms in the Einstein frame; scalar tensor theory 
with matter terms may yield the same results 0, Il9| . 
An interaction term of the form Q = 3cH(p p + p m ) was 
proposed to look for a dynamical solution to the coinci- 
dence problem pp| . Such the interaction term might be 
motivated by analogy with dissipation. For example, a 
fluid with bulk viscosity may give rise to a term of this 
form in the conservation equation. Without an interac- 
tion term, it was shown that the dark matter could not 
track the phantom energy and would be quickly driven 
to zero [jjj • It is of interest to study the cosmological 
consequences of the above two types of interactions in 
the phantom energy model. 



III. MODEL I 



Let us first consider the following interaction [Tsl TTflL 



-CKp m 0, 



(13) 



where c is a dimensionless parameter. Such a coupling 
arises for instance in string theory [Isj . or after a con- 
formal transformation of Brans-Dicke theory pJJ]. The 
evolution equations I|8I1U|I can then be written as an au- 
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tonomous system: 

x 1 = -3x (l + x 2 



V6 



c z 



v/6 



y' = -3y ( x 2 + -^Ax - ^z 2 



-3z 



1 



6 C 



Ay 2 , (14) 
(15) 

(16) 



which has three critical points as follows. 
Point A: 



xa 



V6 



VA 



0. 



ZA = A/1 



This solution is physically meaningless since Q r 
c ^ 0. 
Pomi B: 



(17) 
> 1 if 



x B 



V6 
6 



A, 



A 2 
6 ; 



z B = 0. (18) 



This critical point corresponds to the phantom- 
dominated solution ftp — 1, which always exists for any A 
and c. The effective equation of state, u> c ff = — 1— A 2 /3, 
depends on the slope of the potential. To find out un- 
der what condition this fixed point is a stable solution, 
we study the behavior of small deviations from the so- 
lution. The linearization of system l|14|) - (|16l) about this 
fixed point yields two eigenvalues mi = — (3 + A 2 /2) and 
iri2 = — (A 2 + 3 — cA). Thus the phantom-dominated 
solution is stable for c < A + 3/ A. 
Point C: 



xc 



4 



V6(A - c) ' 



2/c = l- 



A 



A-c (A-c) 2 



2(A-c) 2 ' 

(19) 



This fixed point corresponds to the phantom-fluid scaling 
solution, which exists for 



c < 



A - VA 2 + 12 A + VA 2 + 12 , 3 
2 ° r 2 ^ A+ A? 



Substituting linear perturbations x — > xq + 5x, y — > 
yc + Sy, and z — > zc + Sz about the critical point into the 
system of equations (|*PI|) - (rRjll . to first-order in the per- 
turbations, gives the following two independent evolution 
equations of the linear perturbations: 



8x 



Sy' 



l + 3x 



c 



y 2 c 



V6 



CXr 



Sx 



3xc Vc + -7jr(A- c)y c 



8y, 



-3y c 



A + x c \Sx - iy 2 c 5y. 



(20) 
(21) 
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FIG. 1: The convergence of different initial conditions to the 
attractor solution in the (x, y) phase space for the model I 
with c = —2 and A = 1. 
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FIG. 2: The evolution of the fractional densities of the phan- 
tom field (the dashed curve) and the dark matter (the solid 
curve) for the case in Fig.l. 



The two eigenvalues of the coefficient matrix of the above 
equations determine the stability of the critical point. We 
find that the solution is unstable. 

In the cosmological model with the interaction 1131) 
between phantom field and dark matter, the phantom- 
dominated solution is the only attractor solution in the 
parameter space, c < A + 3/ A. In Figs. 1-4, we plot the 
numerical results. Comparing Fig. 2 to Fig. 4, we see that 
the phantom energy more quickly dominates the universe 
when the parameter c decreases. In Fig.l and Fig. 3, the 
trajectories converge at the same fixed point, which is 
only determined by the parameter A. Hence energy trans- 
fer whether from the phantom field to the dark matter or 
vice versa yields the similar cosmological consequences. 
The stable critical point B with xb < indicates that 
the phantom field climbs up the exponent potential. The 
energy density of the phantom field increases as the uni- 
verse expands, which leads to unwanted future singular- 
ity, and therefore the coincidence problem becomes more 
difficult. In the next section, we will investigate a phe- 
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FIG. 3: The convergence of different initial conditions to the 
attractor solution in the (x, y) phase space for the model I. 
We choose c = 1 and A = 1. 



written as an autonomous system: 
x' — — 3x ( 1 + x 



1 2 1 -2 
-z — -ex 



I) 



-3j/ x 2 



^—Xx z 2 

6 2 



-2>z ( - +x 



2 1 2 1-2 

-z — —cz 



which has four critical points. 
Point A: 



1 



V6 



Xy 2 , (23) 
(24) 
(25) 



= -(VT + ic-l), y A = 0, 



(Vl+Ic+i). 



(26) 




This solution is physically meaningless since fl p < 0. 

Points B,C,D: The other three critical points are so- 
lutions of the following set of equations: 



f(x) = c, 

.,2 



V 



x 2 — ^z-Xx + 1, 



2 , V6, 



z 2 = 2.r' ! + — Ax, 
where we have defined a cubic function 



(27) 
(28) 

(29) 
(30) 



FIG. 4: The evolution of the fractional densities of the phan- 
tom field (the dashed curve) and the dark matter (the solid 
curve) for the case in Fig. 3. 



nomenal model, in which the cosmic doomsday is avoided 
and the coincidence problem may be alleviated. 



IV. MODEL II 



Now let us consider the specific interaction [ljj, |20j, |2 



The critical point with xb < 0, labelled by B, exists for 



< c < /( = ). 

There are two critical points with xc.d > 0, one of which 
is physically meaningless, labelled by D. The other point, 
labelled by C, exists for 



< c < mm{/( j= ), /( — = )}. 



V6 



2V6 



We see that the point B corresponds to a climbing-up 
phantom field, while the point C corresponds to a rolling- 
down phantom field. In order to study the stability of 
the two critical points, we obtain the two independent 
evolution equations of the linear perturbations 



Q = 3cH(p p + p m ), 



(22) 



where c is a dimensionless parameter denoting the trans- 
fer strength. This type of interaction has been proposed 
to look for a dynamical solution to the coincidence prob- 
lem in Ref . [2(| . Then the equation system I|8I1U|I can be 



Sx' = -- [l + 3x 2 B C + y 2 B C + — 



Sx 



b B,C 



3xb,c Vb,c + ~^ x Vb,c ) <%) 



(31) 



Sy' = -3y B C \^X + x B C \5x-3y 2 B C Sy. (32) 
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FIG. 5: Stability regions of the (A, c) parameter space for 
the model II. In the regions I, either the climbing-up scaling 
solution or the rolling-down scaling solution is the stable late- 
time attractor. In the region II, the climbing-up solution is 
the stable late-time attractor. The solutions are physically 
meaningless in the region III. 



The corresponding eigenvalues of the coefficient matrix 
of the above equations indicate that the critical points B 
and C are always the late-time stable attractor solutions 
if they exist. 

In the case of the interaction form ilL'L'l) , there exist two 
kinds of stable scaling solutions, the climbing-up scaling 
solution with xb < in Figs. 6-8 and the rolling-down 
scaling solution with xc > in Figs. 9-11. As shown in 
Fig. 7 and Fig. 10, the universe evolves from the matter- 
dominated phase to the scaling solution, which is charac- 
terized by a constant ratio of the energy densities of the 
dark matter and the phantom field. This may provide 
us with a phcnomcnological solution of the coincidence 
problem. These results agree with those in Ref. [l7| . 
The different regions in the (7, c) parameter space lead 
to different qualitative evolutions. In the region II of th 
parameter space in Fig. 5, the critical point B is a sta- 
ble solution. However, in the region I both the points 
B and C are stable. Which one is the late-time sta- 
ble attractor solution? The phantom field either climbs 
up or rolls down the exponent potential, which is deter- 
mined by the initial velocity of the phantom field. If the 
phantom field initially climbs up, the effective equation 
of state w e ff tends to below —1 and realizes a transition 
from w c s > — 1 to w c s < — 1 in Fig. 8. Thus the universe 
ends with a big rip. If the phantom field initially rolls 
down, the effective equation of state w e s tends to above 
— 1 and realizes a transition from Weg < — 1 to Weg > — 1 
in Fig.ll. In this case the cosmic doomsday is avoided 
and the universe accelerates forever. 



V. CONCLUSIONS AND DISCUSSIONS 

We have presented a phase-space analysis of the evolu- 
tion for a spatially flat FRW universe driven by an inter- 



FIG. 6: The convergence of different initial conditions to the 
attractor solution in the (x, y) phase space for the model II 
with c = 0.2 and A = 1. We choose initial conditions with 
x < 0. 
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FIG. 7: The evolution of the fractional densities of the phan- 
tom field (the dashed curve) and the dark matter (the solid 
curve) for the case in Fig. 6. 




FIG. 8: The evolution of the effective equation of state for 
the case in Fig. 6. 
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FIG. 9: The convergence of different initial conditions to the 
attractor solution in the (x, y) phase space for the model II 
with c = 0.2 and A = 1. We choose initial conditions with 
x > 0. 
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FIG. 10: The evolution of the fractional densities of the phan- 
tom field (the dashed curve) and the dark matter (the solid 
curve) for the case in Fig. 9. 
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FIG. 11: The evolution of the effective equation of state for 
the case in Fig. 9. 



acting mixture of dark matter and phantom field with an 
exponent potential. We have discussed the existence and 
stability of the cosmological scaling solution for two types 
of interactions, namely Q = —CKp m <j) motivated by the 
conformal relationships between the Jordan and Einstein 
frame in string theory, and Q = 3 c H(p p + p m ) motivated 
by analogy with dissipation. In the former model, the 
phantom-dominated solution is the only attractor solu- 
tion when c < A + 3/ A. Energy transfer whether from the 
phantom field to the dark matter (i.e. c > 0) or vice versa 
(i.e. c < 0) leads to similar behavior, except that the 
phantom energy more quickly dominates the universe in 
the latter case than the former case. Since the phantom 
field climb up the exponent potential, the energy density 
of the phantom field quickly increases as the universe ex- 
pands, which leads to an unwanted future singularity. 

However, in the model II with the interaction 12211 be- 
tween phantom field and dark matter, there exist two 
kinds of stable scaling solutions, the climbing- up scal- 
ing solution and the rolling-down scaling solution. The 
existence of a stable scaling solution requires a transfer 
of energy from the phantom field to the dark matter. In 
this model the universe evolves from a matter-dominated 
phase to a scaling solution, which is characterized by a 
constant ratio of the energy densities of the dark matter 
and the phantom field. This may provide us with a phe- 
nomenological solution of the coincidence problem. Fur- 
thermore, in the climbing-up case, the effective equation 
of state w e s tends to below — 1, and then the universe 
ends with a big rip. The effective equation of state may 
realize a transition from w c g > — 1 to w c g < — 1. In 
the rolling-down case, the effective equation of state w e s 
tends to above —1 and may cross —1. In this case the 
cosmic doomsday is avoided and the universe accelerates 
forever with a power-law form. What is the ultimate fate 
of the universe? The universe either accelerates forever 
or ends with a big rip, which is determined by not only 
the model parameters but also the initial velocity of the 
phantom field in the scenario of interacting phantom en- 
ergy with dark matter. 
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